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Abstract 

A new classification scheme for real numbers is given, motivated by ideas from 
statistical mechanics in general and work of Knauf 17 and Fiala and Kleban [9] in 
particular. Critical for this classification of a real number will be the Diophantine 
properties of its continued fraction expansion. 

1 Introduction 

Though this paper is about number theory in general and about a classifica- 
tion scheme for real numbers in particular, it has its roots in Thermodynamic 
FormaUsm, which was developed in the 1960s by Ruelle [32] [33] , Sinai [35] and 
others in an attempt to put statistical mechanics on a firm mathematical foun- 
dation. Once done, the underlying mathematical scheme can then, in principle, 
be applied to non-physical situations, using the original real- world interpreta- 
tions to guide and influence what questions are to be asked and what structure 
is to be discovered. 

This process has been begun in number theory. In [17j . Knauf developed a 
one-dimensional thermodynamic system based on the Farey fractions that ex- 
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hibited phase transition. In [3] , Fiala and Kleban generalized Knauf 's work and 
showed that their generahzation has the same free energy as Knauf 's. We will 
put these earlier works into a common linear algebra framework, allowing us to 
make a seemingly minor, but actually significant, change in the original parti- 
tion function. We will produce, for each positive real number, a thermodynamic 
system. Different real numbers will exhibit different free energies, giving us a 
new classification scheme for positive real numbers. (This classification scheme 
can easily be extended to also include negative reals.) 

In Section 2, we give a brief overview of the parts of the statistical mechanics 
formalism that we will be using. In particular, we will see the key importance 
of the partition function. In Section 3, we tic this formalism to number theory, 
in particular to the Farey matrices. In Section 3.3, we put Knauf's work into 
this language and do the same thing in Section 3.4 for Fiala and Klcban's work. 
In Section 3.5, we show how to alter the earlier partition functions that will put 
us into the world of Diophantine analysis. This is the section in which we are 
not just changing the notation from earlier work. 

In Section 4.1, we use our Diophantine partition function to give a new 
classification scheme for real numbers. In particular we develop the idea of a 
real number having a 1-free energy limit. In Section 4.2, we show how this is 
naturally linked to continued fractions. The rest of Section 4 deals with proving 
that there are real numbers with 1-free energy limits, that there are reals without 
a 1-free energy limit, that all algebraic numbers have a fc-free energy limits with 
k > 1 and that all quadratic irrationals have 1-free energy limits. We also show 
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that e has a \/iV log A''-free energy hmit. We will conclude with open questions 
in Section 5. 

There has been a lot of other work linking statistical mechanics to number 
theory. There is other work of Knauf [T^ [20] [21] , of Guerra and Knauf [12] , of 
Contucci and Knauf [5], of Fiala, Kleban and Ozliik [lOj, of Kleban and Ozliik 
[IB] , of Prellbcrg, Fiala and Kleban ^5], of Feigenbaum, Procaccia and Tel [5] 
and others. 

There is also the transfer operator method, applied primarily to the Gauss 
map, which allows, in a natural way, tools from functional analysis to be used. 
We believe this was pioneered by Mayer (see his j25j for a survey) , and non- 
trivially extended by Prellberg [271, by Prellberg and Slawny [29 , by Isola [2] 
and recently by Esposti, Isola and Knauf 7 . An introduction to this work is in 
chapter nine of Hensley [13]. We will not be following this approach here. 

As of August 2008, the web site 
http://www.secamlocal.ex.ac.uk/people/staff/mrwatkin/zeta/physics.htm 
offers many other attempts over the years to find links between statistical me- 
chanics and number theory. 

Finally, I would like to thank Edward Burger for many interesting conver- 
sations about this work and Steven Miller and L. Pedersen for comments on an 
earlier draft. Also, I would like to thank Peter Kleban, Ali Ozliik and Thomas 
Prellberg for finding a significant error in an earlier draft. 
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2 The Partition Function and the Free Energy 

This is a rapid fire overview of basic terms in statistical mechanics. For each 
A'' e N, we have a finite set Sjv, called the state space. Let 

£ : Sjv ^ M+ 

be a function that we call energy. The partition function is defined to be 

If we were modeling a physical system, the elements in the state space corre- 
spond to what can happen. The variable /3 corresponds to the inverse of the 
temperature. The underlying physical assumption is that the probability that 
a system is in a state a G §n will be 

g-/3E(a) 

Probability in state a = . 

Zn{p) 

While far from a proof, this interpretation makes sense, in that at high temper- 
atures (meaning for (3 close to zero), all states become increasingly likely, while 
at low temperatures the most likely state increasingly becomes the state with 
the lowest energy. 

There is a free energy if the following limit exists: 

N^oo iV 

with the function /(/?) being called, naturally enough, the free energy. It is 
believed that phase transitions occur at values of /3 for which /(/?) fails to be 
analytic. 
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For almost all of this paper, our state space will be 



Sat = {o- = (cti, • ■ • ,0-Ar) : <Ti=0 or 1}. 



Thus each of our S^r will have order 2". We can think of our state space as 
having N site points, each having value or 1. 

The most famous example is the one-dimensional Ising model. For conve- 
nience, we let each site have the value of 1 or —1. Thus for the Ising model, we 
have 



Ising, in his 1925 thesis, showed that for this model there is no phase transition, 
meaning he showed that the free energy is an analytic function. For the two- 
dimensional analog, it is one of the great discoveries (originally by Onsager in 
1944) that phase transition does occur. Most texts on statistical mechanics, 
such as [37] J describe the Ising model in detail. 

Note that in the Ising model, a site will only interact with those other 
sites that are immediately adjacent to it. This is an example of finite range 
interaction. Since there is no phase transition for the one-dimensional Ising 
model, it was long believed that there would be no phase transition for any one- 
dimensional system. But in the 1960s, it was discovered that phase transition 
can occur if the interactions are not of finite range but over possibly arbitrarily 
long distances. A good introduction to this work is in Mayer's The Ruelle-Araki 



Sat = {(T = (di, . . . , (Tn) ■■ <7i = ±1}. 



The energy function for the Ising model is 



N 
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transfer operator in classical statistical mechanics [24j . Such interactions are 
called long range interactions. In the following number theoretic models, it is 
key that the interactions are long range. 

3 Number Theoretic Partition Functions 
3.1 General Set-up 

Fix a positive integer k. For each positive integer N, our state space will be 

Sat = {(ci, . . . , (Tjv) : CTi 0, 1, . . . , - 1}. 

Thus S^v contains elements. 

We define a new type of product of an TV-tuple of rt x n matrices with an 
M-tuple of such matrices to be the MA^-tuple: 

{Ai,...,An){Bu...Bm) = {AiBi,AiB2,...,AnBm). 

For matrices A = {uij) and B = (bij), we denote the Hilbert-Schmidt prod- 
uct (which is also called the Hadamard product) as 

A*B = Tr{AB'^) = ^ aijhj. 

l<i,j<n 

For example, thinking of a 2 x 2 matrix as an element of M^, then A * B = 
Tr{AB^) is simply the dot product of the two vectors. 

Let M„ denote the space of n x n matrices. For a function 

/ : M„ -> R 

and for two n x n matrices M and A, define 
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following notation as in . For k n x n matrices Ai, Ai, . . . , Ak, define 

k 

fiAimui,. . . ,Ak) = j2f{Mm\A. 

Consider the map 

Z:NxM„xMxM^x r(M„,M) K*, 

where N is the natural numbers, M„ is the space of n x n matrices, r(!M„,M) 
is the space of functions from n x n matrices to the real numbers and M* is the 
extended real numbers, defined by setting 

Z{N, M, /3, (Ao, . . . , Ak-i), f) = /(M)(/3)|(Ao, . . . , A^^.f . 

Here the notation (Aq, . . . , Ak-i)^ is referring to the above newly defined prod- 
uct of tuples of matrices and hence can be viewed as short-hand for all products 
A,^ • •• with < < fc - 1. 

We want to link this with partition functions. Fix an n x n matrix M and 
also k n X n matrices Aq,Ai, . . . ,Ak^i. Let the function / be the constant 
fimction 1, or in words, let f{B) = 1(B) = 1 for all matrices B e M„. Then 
define the partition function to be 

Zjv(M)(/3) = ZiN, M, f3, (Ao, . . . , Ak-i), 1) = l(M)(/3)|(Ao, . . . , Ak-i)"" . 

The "physical" intuition is as follows. Think of a one-dimensional lattice 
with n sites. At each site, there are k possible states, each of which can either 
be indexed by a number ai between and fc — 1 or by a matrix A^. . Then the 
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states can be viewed as either all possible 

cr = (cTi , . . . , (T„ ) G Z/fcZ 

or all matrices of the form 

-Aai • • • ^(T„ ■ 

In order to get the above partition function, we set the energy of a state to be: 

E{a)^log\M^{A,,---A,J\. 

In our applications, it is more natural not to emphasize the energy function. 

Using this language, we have natural recursion relations linking the partition 
function Z„_|_i with partition functions for various Z„, with different choices for 
the matrix A'l. More precisely 

Lemma 3.1. 

ZN+i{M){f3) = Zn{MAJ){i3) + ... + Zn{MAI){(3). 

This is a simple calculation. It offers a more general and natural form for 
the key recursion relation (2) in [9J. 

3.2 Farey Matrices 

This section continues the building of needed machinery. (See, also, section 
4.5 in [TT]. Another source would be [53]) 

We develop the Farey partitioning of the extended real numbers. Start with 
the set 
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and define the Farey sum of two fractions, in lowest terms, to be 

p r p + r 
- © - = . 

q s q + s 

We now extend a given nth Farey set S'„ to the {n + l)st Farey set by adding 
to 3^n all of the terms obtained by applying the Farcy sum. Tlius wc have, 

~ to'T'T 

1 2 1 1 01 

o'T'T' 2' ij 

_ fl 3 2 3 1 2 1 1 01 
^ ~ to'T'T'2' 1'3' 2'3'Tj ■ 

By reversing the order of the terms in we get a partitioning of U oo. 
Here we are thinking of ^ as the point at infinity, which is why we are working 
with the extended real numbers M* . 

We now describe this partitioning in terms of iterations of matrix multipli- 
cation. Let 

^o=(; J) and^i=( J I). 

These two matrices are key to understanding the Farey decomposition of the 
unit interval and continued fraction expansions. Further, these two matrices 
will be key to the partition function of Knauf [17] , of Fiala and Kleban <^ , and 
to our eventual use of partition functions for Diophantine approximations. 
Note that 

f P \ A _ f P + T ^ ^ 
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P r \ . I p p + r 
A\ = 



and 

( 9 s y ' V « ^ + •5 

Then 3"!, save for the point at infinity ^, can be obtained in the natural way 
by examining the right columns of first then ^o- Likewise, 3^2, save for the 
point at infinity can be obtained in a similar natural way by examining the 
the right columns of first then AiAq, then AqA\ and finally Ag. In general 
the elements of "Jn^ save of course for the point ^, are given by the right columns 
of the 2" products of matrices of A\ and Aq. 

These allow us to recover the continued fraction of a positive real number 
a. We know that any real number a can be written 



1 

a = ao + 



ai 



T 03 + . ■• 

which is usually denoted by 

a= [ao;ai, 02,03,...], 

where ao is an integer and the remaining are positive integers. The number 
a is rational if and only if its continued fraction expansion terminates. We say 
that the rational is the mth partial fraction for the number a if 

Pm r -1 

— = [ao; ai, . . . , (imj. 

We want to use our Farey matrices to find the nth partial sum 2^. for a 
given number a. We return to our Farey numbers, but reverse the orders of the 
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numbers: 

= {t't4} 

T = JO i 1 2 n 

^ ll'2'l'l'Oj 

_ ro 1 1 2 1 3 2 3 n 

^ ~ \1'3'2'3' 1'2' 1' I'Oj" 

We can thus view J'l as providing a splitting of the positive reals into two 
intervals: and [j,^]- Note that can be thought as flipping the 

columns of Aq and [j, can be thought of as the flipping of the columns of A-i. 

Likewise, 9^2 will split the positive reals into four intervals: [j,-^] (which 
can be thought of as the flipping of the columns oi AqAq), (which can be 

thought of as the flipping of the columns of AqAi), [j, j] (which can be thought 
of as the flipping of the columns of AiAq) and [j, ^] (which can be thought of 
as the flipping of the columns of AiAi). 

This pattern continues. Consider the matrix 

Aao Aai aO-n-i aun 

with ao > and o, > 0. Then the left column of the matrix will correspond to 
[ao; ai, . . . , ajv] while the right column will correspond to [ao; ai, . . . , cijv-i]- 
For matrices 

Aao /jai ^ojv-i AaN 
^1 ^0 ■ ■ ■ ^0 ^1 

with ao > and a, > , then the left column will correspond to [ao ; a i , . . . , cat _ i ] 
while now the right column will [ao; ai, . . . , ajv]- 
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Thus to determine the continued fraction expansion for a given positive real 
number a, we just have to keep track of which interval a is in for a given Jjq. 

3.3 Knauf's work 

This section will show how Knauf's number theoretic partition function [17j 
can easily be put into the language of this paper. 
Consider the sets J'jy H [0, 1] 



As is well known, as — + cxd, these sets will eventually contain all rational 
numbers in the unit interval [0, 1]. Then Knauf defines his partition function as 



f e;rNn(0,i) 

(Note that the K is not being used as an index but stands for 'Knauf; the 
subscript N is an index.) In the limit we get 



where ^(n) is the Euler totient function. In turn, X^J^i <l>(!n)l'n^ is well-known 
to equal to 






n=l 



C(/3-l) 

m 
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for (3 > 2, where CiP) is the Riemann zeta function and is not defined for P <2, 
showing that there is critical point phenomena for this one-dimensional system. 
The C(/3 — l)/C(/5) '^ill show up throughout this paper and is why many of the 
later theorems are only true for f3 > 2. 

We now show how this can be put into the language of Section 4.1. Let 




1/' 



and let Aq and Ai be the above Farey matrices. 
Then the Knauf partition function is 

Z^(/3) = l(M^)|Ao(Ao,Ai) 



N 



The initial Aq is just to insure that we are in the unit interval. Also, for 
Ao{Ao, Ai)^ , we are using the new product for matrices defined in section 3.1 
and not traditional matrix multiplication. 

3.4 Fiala-Kleban Work 

In [5], Fiala and Kleban considered a different number theoretic partition 
function. In the language of this paper, let 



M 




X 1 



(Again, the F is not an index but stands for "Fiala-Kleban".) Then the new 
partition function will be 



Z^{x,P) = \{M^)\Ao{Ao,A^) 



N 



As mentioned earlier, the recursion relation (2) in Section 2 of [9] is simply a 
special case of Lemma 3.1. 
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Two partition functions are said to have the same thermodynamics if their 
free energies are equal. Thus, the Knauf partition function Z^(/3) and the 
Fiala-Kleban partition function (/3) have the same thermodynamics if 

li^ log(^f(/?)) ^ log(^^(/3)) 

N^oo N N~*oo N 

This equahty is shown in Section 4 of [9J , using as an intermediary tool a certain 
transfer operator and depending on the earlier work of Knauf |17j . In [lOj , Fiala, 
Kleban and Ozliik showed the thermodynamics of the Knauf partition function is 
thermodynamically equivalent to a number of other number theoretic partition 
functions. It is certainly the case, though, that thinking of the various matrices 
as vectors in will yield more straightforward proofs of these equivalences. 

3.5 A Diophantine approach 

We now make a seemingly minor change in our choice for the matrix M that 
will create a quite different thermodynamics, leading in the next section to a 
new classification of the real numbers. Set 



M 



-1 
a 



for some real number a. Define the Diophantine partition function to be 

ZN{a;l3) = l(M)|(Ao,Ai)^. 



Then we have 
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Note that the terms that dominate the above sum occur when |p — agj is smalL 
This places us firmly in the realm of Diophantine analysis. We will see that we 
can classify real numbers a by understanding the existence of free energy limits 
for the statistical systems associated to the partition function Z^la; (3). Again, 
while this partition function is cast in the same overall language as Knauf and 
Fiala and Kleban, its thermodynamic properties will be quite different. This is 
what separates the present work from [17] and [9]. 

4 Classifying real numbers via free energy 
4.1 The classification 

With the notation above, we have for each positive real number a and each 
positive integer N the partition function Ziq{a; (3). 

Definition 4.1. A real number a has a k-free energy limit if there is a number 
Pc such that 



exists for all (3 > fic 

For fc = 1, this is a number-theoretic version of the free energy of the system. 
By an abuse of notation we will also say that a has a /(iV)-free energy limit, 
for an increasing function f{N) if 



exists for all j3 > fic- 

To see that this is a meaningful classification scheme for real numbers, we 
will establish the following three theorems: 



lim 

Af— >oo 



log(Zjv(a;/3)) 



lim 

N^oo 



log(Zjv(a;/3)) 
fiN) 
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Theorem 4.2. There exists a number that has a 1-free energy limit for /? > 2. 

Theorem 4.3. There exists a number that does not have a 1-free energy limit, 
for any value of (3. 

Theorem 4.4. Let a be a positive real number such that there is a positive 
constant C and constant d>2 with 

C , 

for all relatively prime integers p and q. Then a has a k-free energy limit for 
/3 > 2, for any k > 1, and in fact, the k-free energy limit is zero. 

An easy consequence of the above is that aU algebraic numbers have k-free 
energy limits equal to zero, for any k > 2. 

While we do not know if algebraic numbers have 1-free energy limits, we will 
show 

Theorem 4.5. All quadratic irrationals have 1-free energy limits for (3 > 2. 

We will also show, in Section 4.7 that the number e has \/N log{N)-iiee 
energy limit. 

4.2 Links to continued fractions 

The goal of this paper is not just to give a new way to classify real numbers 
but also to show how such a classification scheme follows from the thermody- 
namical formalism, fitting into a more general framework. But if all we wanted 
was the classification scheme, then it is possible to reframe our definitions so 
that there is no need for the language of statistical mechanics. The goal of this 
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section is to state the theorems that would allow us to avoid thermodynamics. 
They will also be key to proving the five theorems of the previous section. 

Let our positive real number a have continued fraction expansion 
[ao; ai, a2, • • •]• We know that the best rational approximations to a are given 
by the rational numbers [ao; ai, 02, . . . , a^,]. The fractions [oq; ai, 02, fc], 
with < A; < 0^+1 , are called the secondary convergents to a. For a given m, 
we know that all the vectors corresponding to the [ao; ai, 02, . . . , am-i,k], with 
< k < am, are on the same side of the line x = ay, while [ao; ai, 02, . . . , am] 
jumps to the other side of the line. 

We now set some notation for the rest of the paper. Given the positive real 
number a, for each positive integer N there is associated a positive integer m 
and nonnegative integer k such that 

N = ao + ai-\ \-am + k, 

with < fc < Om+i- We create a subsequence of N, denoted by Nq, Ni, N2, . . ■ 
by setting 

Nm = ao + ai H 1- am- 

In this notation, given any A'' e N, we have 

Nm<N< Nm+U 

or, in other words, 

ao + ai H h am < A'' = ao + ai H h a^ + A; < oo + ai H h a^ + a^+i . 
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Set 

Pn , ,1 

with < fc < a„i+i, with pjv and ^at having no nontrivial common factors. We 
know that the fractions are the best rational approximations to the initial 
real a. Denote 

dN = r- 

\Pn - aqN I 

We have the following chain of inequalities that will be key: 



dN^-i < dNm + l < dN^+2 < ■ ■ ■ < dAr^+a„+i-l < dN^, 

which are well-known (for motiviation, see the chapter on continued fractions 
in [5S]). 

We will show 

Theorem 4.6. For any (3 > 2 and for any positive real number a, we have for 
all positive integers N that 

log(4) ^ log(Z^(a;/3)) ^ log(^iV4„) 
N - N - N 

To show that there are numbers that do not have 1-free energy limits, we 
will construct an a so that a subsequence of approaches infinity. In 

turn, to show that that there are numbers that have 1-free energy limits, we 
will construct an a so that iW)^^'^^,,,) ^ppj-oaches zero, forcing '°s(^J^("-^)) 
to approach zero (and thus guaranteeing that the limit exists). 
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4.3 Preliminaries 

The reason why we look at the products of the matrices Aq and Ai is that 
the right columns will correspond to all of the integer lattice points ( g ) with p 
and q relatively prime in the the first quadrant of the plane. If Ai^ ■ ■ ■ Ai^ has 
right column ( ^ ) , then 

M*{A,,---AiJ = {-l,a)- Q = aq-p 

Let 

be vectors that satisfy dot (wi, 172) = ±1. Let C(wi, W2) denote the cone of integer 
lattice points defined by: 

C{vi,V2) = {avi + bv2 ■ a, b relatively prime nonnegative integers}. 

Let Cjv(fi, ^2) be the subset of C{vi,V2) consisting of vectors that are the right 
columns of all possible Ai^ ■ ■ ■ Aj„ . Then set 

ZN{a;P,Vi,V2) = -. ^— r^. 

In the same way, we set 

Z{a;P,vi,V2) = V -. — TT. 

Suppose we have integer lattice vectors vi = (i) ,V2, ■ ■ ■ ,Vm on one side of 
the line x = ay and integer lattice vectors wi = [q),iu2, ■ ■ ■ ,'Wn on the other 
side of the line x = ay such that det(wj, t^i+i) = —1, det(u'j, Wj+i) = 1 and 
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det{vm,Wn) = — 1. Then we have 



m — 1 



n-1 



(a;/3)< ^ZN{a;f}, 

i=l i=l 

We need the above to be an inequaUty since there is "overcounting" on the 
right hand side, since, for example, the part of Z]\[{a; P,Vi-i,Vi) coming from 
the vector Vi also appears as a term in Z^la; /3, Vi,Vi+i). The key, as we will see, 
is that the Z^ia; /3, Vm,Wn) term contributes the most to the partition function 
Zffia; 13). For the rest of this section, we want to control the sizes of the various 
ZN{a;l3,Vi,Vi+i) and Z^ia; P,Wi,Wi+i). 

We first return to the more general case of two vectors t^i = ( ) , i'2 = ( ) 
with det(wi,W2) = ±1, under the additional assumption that vi and V2 lie on 
the same side of the line x = ay. 

For 



di = 



1 



i,d2 = 



1 



|pi - agil \p2-ccq2\' 
suppose that di < d2, which means that the line through the origin and the 

point {p2, (72) is closer to the line x = ay than the line through the origin and 

We want to show that 

Z{a;p,v^,V2)<\d2\^^^^^^. 



We know that 



1 



(-l,a) 



1 

^2 



Let V be some integer lattice point in the cone C{vi,V2). Then there are 
relatively prime positive integers a and b with v = avi + bv2. 
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We have 

{a + b)\{-l,a) ■V2\ < \{~l,a) ■ {avi+bv2)\ 

< (a + 6)|(-l,a) -wil. 
Inverting and raising everything to the power of f3, we have 

{a + b)\P ^ \{-l,a)-v\P ^ (a + 6)p- 
Summing over every vector in C(wi, 112), we get 

where in the first and third summation we are summing over all relatively prime 
positive integers a and b. It is well-known, as mentioned in Section 4.3 of Knauf 
[TT] . that, 

1 C(/3~l) 
2^\a + b\P C(/3) ■ 

It is here that the ^^^^^^^-^ makes its critical appearance. We have our desired 
inequality. 

4.4 Proof of Theorem 4.6 

First, the partition function Zi^[a;(3) is the sum of many positive terms, 
including dP^. Thus we have the lower bound 

log(4) < log(^jv(«;/3)) 
N - N ' 

Now for the upper bound. We are at the Nth stage. Letting t < N, express- 
ing as = flo + ■ • • + a„j + fc, with < fc < am+i and using the notation 
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from section 4.2, we assume that the vectors (qlzl) and (g') He on the same 
side of the hne (x — ay). This is equivalent to there being some s < m with 
ao + -- -+as<t— l<t<ao + -- - + a^+i- We know that 



det(^*-i Pn^±l. 
qt-i qt 



We have 



qt-ij \qtjj V \qt-iJ \qt 

cm 



using that dt < dAr,„ . 

Since there wiU be N such cones, we have 



iV - TV ' 

finishing the proof. 

4.5 Proof of Theorem 4.2 

We know that the best rational approximations to a real number a are the 
— ao; oi, • ■ • , am • 

It is well known that 

9Ar„+i = am+i9Ar„ + qNm-i- 
Further (as in Lemma 7.2 of ^) 

|PAf„ - agjv,„ I 
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Our goal in this section is to construct a real number a for which 
lim ^ — = 0, 



for P > 2, which by Theorem 4.5 will force a to have a 1-free energy limit. 
If 



log(Ari§zl)(rf^j/5) 



JV->oo N 

exists, then it equals 

C(/3-l) 



lim ^+ lim l^i#r+ lim ^'^^'^^ 



The first two terms in the above certainly go to zero. Thus we must construct 
a real a so that limjv-»oo ^ '"^"^ = o. Recalling our notation that Nm = 
ao H + dm < N = Nm + k, with < fc < am+i, we have for each A'' 



N - Nm 

Thus all we have to do is construct an a so that limjv^oo ^ '"^l""* ~ ^• 
For a = [ao; ai, . . .], define the function f{m) by setting 



_ /(m) 
O-m+l — Qn ■ 
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We have 



/31og(rfjv„) ^ /31og(djv„ 



N - 

(3 log[(o„+i + 2)gjv„ 



< 



< 



/31og(2a„+iq'iv^) 



/31og(2) /?l"s('/A'r^'^^^0 



/31og(2) ^ /3[/(m) + l]log(gjv^) 

Since the first term in the last equation goes to zero as A'' — > oo, we have, if the 
limits exist, that 

PlosjdNj ^ f3[f{m) + l]\og{qNj] 

JV^oo Nm ~ iV->oo Nm 

We now start with a go and a qi and a function /(m) and use these to create 
our number a, Let qo = 1, qi = 2, and /(m) = m for m > 1. Then define for 
m > 2, a„+i = gj;^ ^ 

Now 

< log(2a„gjv„_J 

= log(2<z/(™-i)+i) 

= log(2) + (/(m-l) + l)log(gjv_J 

Then we have 



/3[/(m) + 1] log(g;v^) ^ /3[/(m) + 1] log(2) ^ /3[/(m) + l](/(m - 1) + 1) log(g;v„ 



iVm " Nm Nrr 
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We will now use that 



Nm = ao + . . . + a„ 



> an 



f(m-l) 



Then 



/3[/(//0 + 1] l»g(ryA:,J I3{m + 1) log(2) p{m + l),u l.^gjrM,,., J 

which has limit zero as A'^ ^ oo. Thus with the choice of the function f{m) — m 
we have constructed a real number that has 1-free energy limit, finishing the 
proof of Theorem 4.2. 

4.6 Proof of Theorem 4.3 

We use from Theorem 4.5 that ^2514) < M^^HiM ^e will construct a 
real number a so that 

lim = oo. 

This will mean that limjv^oo will not exist, which is the goal of The- 

orem 4.3. 

Proceeding as in the proof of Theorem 4.2, we define a = [oo; ai, 02, . . .] 
inductively on the a'^s by setting, as before, ao = 1 and Um+i = Qn™^ but now 
defining /(m) as 

/(m) = ao H \-am- 
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We use that am+i9iv„ < c^at^ and = cto H + cim, we have 



/31og(<zj;^'g;v.) 



/3(/(m) + l)loggjv^ 

^ /?(ao H am + l)loggiv„. 

oo H ha™ 

Since the last term goes to infinity as Nm — * oo, we are done. 

4.7 Proof of Theorem 4.4 

We assume that a is a positive real number such that there is a positive 
constant C and a constant d>2 with 

^ < \p-m\ 

for all relatively prime integers p and q. 
In particular, we have 

< a . 



In our notation, this is 



where c is a constant depending on a but not on Nm- In a similar argument 
as in Theorem 4.2, the number a will have a k-free energy limit if we can show 
that 

lin, i#^=0. 
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This will certainly happen if we can show 



and hence if 



We know that 



Then we get 



= amqNm-i + lN,n-2 

= 2am{cim.-l(lN,n-2 + iNm-s) 

< 2'^amam-lQNm-2 

< • • • flo- 

Thus we want to examine 



log(2™+iamam-i • . • ao) 

lini ^ TTT-TZ 



or 



(m + 1) log(2) + log(ao) + . . . + log(a„) 
nm r-r . 

Nm^co {ao + ... + ttmr 

Since it is always the case that uq + . . . + am > m + 1, the above limit must 
always be zero for any A; > 1. 
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Corollary 4.7. All positive algebraic numbers of degree greater than or equal 
to two have k-free energy limits equal to zero, for any k > 1. 

Liouville's Theorem (see Theorem 1.1 in [5]) states that all irrational alge- 
braic numbers a have the property there is a constant d> 2 with 



for all relatively prime integers p and q, allowing us to immediately use the 
above theorem. 

4.8 Quadratic Irrationals have 1-free energy limits 

The key will be that every quadratic irrational has an eventually periodic 
continued fraction expansion (see section 7.6 in [36J. We will show first, for a 
quadratic irrational, that lim,„.^oo log(7Ar„/iV„j exists and then show that the 
existence of this limit is equivalent to the number having a 1-free energy limit. 

Let a be a quadratic irrational. We can write 



— <\p- aq\ 



a = 



[bo . . . ,bp,ci,C2, ■ ■ ■ ,c/,ci,C2, . . . ,c/, . . .]. 



The period length is I. For notational convenience, set 



b 



bo + --- + b., 



'p 



c 



ci H h Q 



B 



C 
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For each positive m> p, there are unique integers A; > and n with <n <l 
such that 

m = p + kl + n. 

For this n, set 

dn = ci H h c„ 

= f 1 y.Yo 1 

\l Cn ) V 1 Ci 

Note that there are only a finite number of possibihties for the dn and the Dn- 
From section 7.6 of Stark, we have 

Now Nm = b + kc + dn grows linearly with respect to k. Since qn^ is the lower 
left term of the matrix DnA'^B, we also have that loggjv^ grows linearly with 
respect to k, with constant leading coefficient. Thus limm^oo must exist. 

Now to show why this will imply that a has a 1-free energy limit. Since 
rfjv^_i < rfjv < rfjv^ and since Nm < N < Nm+i, we have from Theorem 4.6 
that 

l°g(4^-J ^ log(^iv(a;/3)) ^ log(^^4^) 

From our earlier work, we know that a will have a l-£ree energy limit, for /3 > 2, 

if we can show 



logrfAT^.i ,. \ogdN„ 
lim = hm 



We will see that these hmits will exist if cind only if liniyj^^QQ 

\r " exists. 
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Let a = [ao;ai,a2, . . .]. We know that Nm+i = Nm-i + + Om+i- For 

m > p, Qm and Qm+i can only be chosen from finite set {ci, . . . c/}. Thus 

logdAr,„_i logdjv„_i 
nm = hm 

m^oo Nm+1 rn^oo Nm-1 

since the denominators, for aU m, differ by a fixed amount (and of course since 

N,n -* OO). 

We know that Qm+i^N^ < rfjv„, < (om+i + '^)qn^- Then, if the limits exist, 
we have 

logarn+iQN^ ^ logdjv™ ^ ,. log(a,„+i + 2)gAr^ 
lim < hni < lim . 

Since the flm+i are bounded and the N.^ — > cxd, we have 

log«,»+i ^ ^.^ log(a„.+i + 2) ^ 

giving us that 

logdAr,„ ,. loggjv™ 
lim = iim , 

concluding the proof. 

Corollary 4.8. The golden ratio <f) = has l-free energy limit equal to 

\og(j), for (i> 2. 

For background on 0, see Chapter 1.7 in [30]. The key is that continued 
fraction expansion for (j) is 

</)=[!: 1,1,1,1,...]. 

From the above theorem, we know that the l-free energy limit, for /? > 2 is 

loggw™ 
lim . 



4 Classifying real numbers via free energy 



31 



For (j), we know that = m + 1 and that qn^ = F^, where is the mth 
Fibonacci number. We know that 



1 / / l + \/5\ / 1-V5) 



V5 U 2 
Then the 1-free energy hmit will be 



lin, 1^ = hm 



m— >oo m— »oo m + 1 

= log?!', 

as desired. 

Of course, similar arguments can be used to find the l-free energy limits for 
many quadratic irrationals. 

4.9 The number e has V^log(Ar)-free energy limit 

We will actually show this for e — 1. The continued fraction expansion for 
e — 1 is 

[1; 1,2, 1,1,4,1, 1,6,...]. 

Thus for e — 1, we have 

ask = ask+i = 1 and 03/5+2 = 2(fc + 1). 

Let qn^ denote the denominators of the fraction [ao;ai, . . .am] associated to 

e - 1, where = ao H H am- 

In an analogous fashion to the proof that quadratic irrationals have 1-free 
energy limits, we will first show for e — 1 that 
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exists and then show that the existence of this Hmit will give us our result. 
We know that 

Then we have 

aoai • • • Um < QN„t < '^"^^^O'mlm-l ■ ■ ■ ClO- 

Then, provided the limits exist, we have 



logaoai • • - an, 
lim ,__ — . < lim 



log qN„ 



< lim 



log2™+^amam_i ...ao 



JVm-cx, VA?ml0g(iV™) JV^-oo V^l0g(iV„) JV^^oo VA^log(A^„) 

We will show that 

logaoai- ■ -am ,. \og2™+^ amam-i ■■ -ao 

hm , ; = hm r= ; ; , 

JV^^oo /?^log(7V„) N^^oo ^/T^log{Nm) 

which will force limjv^^oo log^jv^/V-^m log(iVm) to exist. 

We now have to look at the explicit values for the am- There are three cases, 
depending on if m is 0, 1 or 2 mod 3. We will let m = 3A; + 2 and show that the 
above limits exist and are equal. The other two cases are similiar. The above 
denominator is 

1 + 1 + 2 + 1 + 1 + 4 + 1 + . . . + 2(fc + 1) = 2(fc + 1) + (fc + l){k + 2), 

which has quadratic growth in k. The numerator on the left hand side of the 
above limit is 

logaoH loga^ = log 1 + log 1 + log2 + log 1 + log 1 + log4 H log2(A: + l) 

= A; log 2 + log 1+ log 2 H log(A; + l), 
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which has fclogfc growth. Similarly for the right hand side of the above limit, 
we have 

log 2™+' 

o-mdm-i ■ ■ -ao = log(aTO+i + 2) + (m+l)log(2) + log(ao) + - • ■+^og{am), 

which also has A; log A; growth, with the same leading coefficient. The above 
limits will exist. 

Now to show that 

^.^ log(ZAr(a;/3)) 
N^oo y/N log TV 

exists. 

In an argument similar to the proof of Theorem 4.6, we have 

log(<„_J ^ log(Z^(a;/3)) ^ log(^^4„) 



^JNm+l\ogNm+\ \/NlogN //V;;iogAr„ 
Following the proof in the last section, we see that 



lim . ^"g^^-- lim 



since N^+i — Nm-i grows at most linearly with respect to m and since each 
Nm grows quadratically with respect to m. 

Since am+iQNm ^ ^Jv^ ^ {o-m+i + '^)QN^(lNmi if t^^^ limits exist, we have 
j.^ log am+igjv„ ^ log ^ j.^ log(a„+i + 2)qN^ 



^Nm log Nm rn^<^ V^m log Nm rn^'x V Nm log Nm 

Since the Qm+i grow at most linearly with respect to m and the Nm grow 
quadratically, we have 

log log(am+i + 2) 
lim — = nm — r= = L), 

m^oc y/Nm log Nm V -^m log Nm 
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giving us that 

lim —= = lim —= , 

log Nm log iVm 

concluding the proof. 

5 Conclusion 



There are a number of directions for future work. It would certainly be 
interesting to know if there are algebraic numbers besides the quadratics that 
have 1-free energy limits. Also, what type of free energies are there for various 
other real numbers. We have shown that some real numbers have 1-free energy 
limits while others do not. This too is just a beginning. What types of limits 
are possible? Are there numbers a for which the sequence '"s^-^^^"'^)) has any 
possible limit behavior? For example, certainly we should be able to rig a so 
that any number can be the limit of the sequence. In fact, we should be able to 
find such sequences with two accumulation points, three accumulation points, 
etc. All of these should provide information about the initial real number a. 

Once we know that there is a free energy limit, then the most pressing 
question is to find for which values of (3 is the free-energy non-analytic. These 
points will be the analogs to critical point phenomena in physics. For the Knauf 
approach, there have been a number of good papers (such as in [17] [18] [19] 
[2Q][2T][2][6][2][10][12][I1][15][22]) exploring the nature of these critical 
points. For our Diophantine partition functions, these types of questions seem 
to be equally subtle. 

The various partition functions Zn{P) depend on the choice of the two-by- 
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two matrix M. Different choices of M give rise to different thermodynamics. 
Knauf's and Fiala's and Kleban's are of the form 



But certainly other M can be chosen and studied. 

Also, Zm{P) depends on the choice of the matrices and Ai. Other choices 
for these matrices will lead to different thermodynamical systems, each with its 
own number-theoretic implications. 

Of course, why stick to two-by-two matrices. This leads to multi-dimensional 
continued fractions. There are many different multi-dimensional continued frac- 
tion algorithms. See Schweiger [34] for a sampling of some. Also, Major [23 has 
some preliminary work on this. Each of these will give rise to a thermodynamical 
system, again with meaning in number theory. 

Finally, there is the question of putting these results into the language of 
transfer operators. (See [3T], [1] for general references.) 
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